Abstract. The frequency-period distribution and other 
properties of Cepheids in IC 1613 are discussed and com- 
pared with those of stars in our Galaxy (Milky Way), 
LMC, SMC, M31 and M33. Taking into account the obser- 
vational limitations and related incompleteness, it is con- 
cluded that the frequency-period distribution of Cepheids 
in IC 1613 is similar to that of SMC; we suspect that a 
much larger number of stars exist in IC 1613 with a period 
of less than 2 d that have not yet been detected. A dis- 
■ \ ■ cussion of the deficiency of fundamental mode Cepheids 
. with periods in the range 8 - 10 d in the Milky Way, M31 
' and M33 is reported. The present data are not sufficient 
CN| . to verify if this is produced by a real bimodal frequency- 
^ ' period distribution or whether depends on the lack of pul- 
Q-i. sating stars in such a period range due to pulsational sta- 
' bility reasons. Some arguments are presented in favor of 
I . a bimodal distribution that is a function of the average 
^ ; metallicity. The Milky Way, M31 and M33 have the two 
, maxima located at the same periods, about 5 and 13 d, 
respectively. A comment on very long period Cepheids is 
also given. 
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1. Introduction 

Cepheids are variable stars that are used to measure dis- 
tances of galaxies in the Local Group and nearby clusters, 
and are the primary calibrators for the secondary standard 
candles that are applied at much greater distances. How- 
ever, they are important also for testing the theories on 
the internal constitution of stars and stellar evolution. In 
particular, the study of Cepheids in nearby galaxies is im- 
portant to understand the effects of different metallicities 
and corresponding mass-^luminosity relations on the pul- 
sational characteristics. The purpose of our project was to 
obtain good light curves of Cepheids in IC 1613 and NGC 
6822. In order to make best use of telescope time and reach 
fainter magnitudes, our strategy was to observe in white 
light, i.e. without filter {Wh photometry). The results of 
the survey for variable stars in four fields of IC 1613 were 
reported in the previous papers of this series (see Man- 
tegazza et al. 
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Paper IV, and references therein), 
while the properties of population I and II Cepheids ob- 
served in a field were discussed by Antonello et al. (1999b, 
Paper II). Recently, the first results for NGC 6822 have 



also been presented (Antonello et al. 2001, Paper V). In 
another paper (Antonello et al. 2000b), the light curves 
of the Cepheids with period near 10 d of IC 1613 were 
analyzed and compared with those of the corresponding 
objects in the Galaxy and in the Magellanic Clouds. In 
the present paper we summarize the characteristics of the 
Cepheids observed in IC1613 and make a comparison with 
other galaxies containing a large number of such stars, tak- 
ing into account the results of other recent surveys. We will 
update the discussion about observational evidence for a 
dip in the frequency-period distribution near 8 - 10 d, 
that is seen preferentially in metal rich galaxies. Presently 
there are three different interpretations of such dips. 

1) In their extensive study on the frequency-period dis- 
tribution of Cepheid variables, Becker et al. (1977; here- 



after BIT) concluded that it is not compatible with a stan- 
dard birthrate function, and can only be explained if an 



ad hoc two-component birthrate function is adopted. The 
primary component is a time-averaged "background" rate 
which reproduces the main peak of the distribution (e.g. 
that in the Milky Way). The second component is pro- 
duced by recent star formation (e.g. in OB associations) 
and is characterized by some cutoff below a critical mass; 
this component reproduces the secondary peak and the 
cutoff explains the dip. 

2) Chiosi ( 1989| ) noted that, since the blue band of 
core He-burning of models with relatively large mass 
(^ 9Mq) tend to move back to the red giant region (a 
tendency enhanced by increasing metallicity), such masses 
would spend a large fraction of the blue band lifetime 
within the instability strip; therefore the excess of long P 
Cepheids would be intrinsic to the stellar models and a 
two-component birthrate is not needed. 

3) Buchler et al. (1997) suggested that the dip could 
be a real deficiency of stars caused by the instability of 
the nonlinear fundamental pulsation cycle in this period 
range, and such stars that cannot pulsate in the funda- 
mental mode actually pulsate in the first overtone one. 

We recall a wish expressed more than twenty years ago 
by BIT (p. 649) in relation to the Cepheids in the seven 
galaxies. Milky Way, LMC, SMC, M31, M33, IC 1613 and 
NGC 6822: "The Cepheids in M33, NGC 6822 and IC 1613 
have only begun to be studied, and much more work needs 
to be done. Further information on the frequency-period 
distribution of any of these seven galaxies would be very 
helpful" . 

2. Cepheids in IC 1613 

2.1. Generalities 

Details of the Wh CCD observations (performed 
with the ESQ Dutch 0.9 m telescope) and reduction 
methods were reported in the previous papers of this 
series, along with the lists of variable stars and their 
coordinatae (Antonello et al |1999a|; Antonello et al. 
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pOOOa ; Mantegazza et al. pOOl ). The finding charts of 
the variable stars are available upon request to the 
authors, or can be found in the Archive at the site 
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Fig. 1. Lower panel: PL relation for IC 1613 Cepheids. 
Crosses: fundamental mode; filled circles: first overtone 
mode; open squares: uncertain mode; open triangle: sec- 
ond overtone candidate; cross in circle: fundamental mode 
Cepheid with a probable companion. For a comparison 
with other surveys, we remark that V ^ Wh + Q.3 for stars 
with cepheid color. Upper panel: full amplitude against P 
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The analysis of the Wh light curves was performed as in 
Paper II, and the Fourier parameters were used to identify 
the pulsation mode of shorter P Cepheids. The mode 
classification of the objects for each field was reported 
in the previous papers; we have however revised it by 
using homogeneous criteria, and in particular we tried 
to classify the stars previously reported with uncertain 
modes. The result is that the F mode pulsators include 
also the following stars: V0236A, V2991B, V3277B, 
V3782B, V0793D, V1307D, and the lOT pulsators in- 



clude also the following stars: V0178A, V0524A, V1296A, 
V1289A, V0223B, V0987B. We recall that, owing to the 
uncertainties on the parameters, the mode identification 
for the shorter P stars is only preliminary. 

In the lower panel of Fig. |l| the PL relation is shown 
for the whole set of stars, which includes 82 F mode, 34 
lOT mode, 3 suspected second overtone mode pulsators, 
and 5 Cepheids with uncertain mode. In the upper panel 
we show the diagram of amplitude vs. P. The three F 
mode Cepheids marked with a cross in a circle have the 
lowest amplitude and are the brightest stars among those 
with similar P; we suspect that they are binary systems 
or there is a blending effect. An analysis of all the Fourier 
parameters indicates that, for P less than about 4 d, the 
formal errors can be relatively large, and the resulting re- 
lation between the phase difference <j>2i = 02 — 20i (of the 
second and first Fourier component) and P is relatively 
scattered. Therefore, this does not allow us to compare 
reliably such parameters of short P fundamental mode 
Cepheids in IC 1613 with those of similar stars in other 
galaxies, when looking for slight differences due to differ- 
ent metallicities. In other words, with the present data, 
useful comparisons shou ld be l imited to stars with P <; 4 
d (e.g. Antonello et al. 2000b ), instead of attempting a 



general comparison, as done in Paper II. 



2.2. Comparison with OCLE survey 

Recently, the OGLE project made a survey of Cepheids 
in IC 1613, in the V and / bands (Udalski et al. ^OOlD , 
using a 1.3 m telescope. Owing to the larger field of view 
(14f), they observed almost the whole galaxy and found 
about 135 population I Cepheids, with a limiting average 
magnitude V ~ 22.6. In the V band the total exposure 
time was 900 sec, and the seeing conditions were much 
better than those of our project; in fact the Dutch 0.9 m 
telescope suffered from focusing problems. In our four se- 
lected fields we found about 128 Cepheids, with a limiting 
magnitude of ^ 23. We verified that, in the same four 
fields, OGLE detected about 75 Cepheids; there are 64 
Cepheids in common with us. The comparison shows that 
there are systematic differences of some arcsec in the dec- 
lination and right ascension between OGLE and our lists. 
The P values are reasonably in agreement. The largest 
discrepancies were found for three stars with P close to 
2 d, in which case it is difficult to distinguish between 
the true P and the alias produced by the one-day sep- 
aration of observations. Indeed, in the power spectrum 
there is a significant peak at both / = 1/P and 
/ = 1 - 1/P d^^. We merged our and OGLE dataset 
after correcting for the different mean magnitude and dif- 
ferent amphtude, and we got for V0122B=OGLE120U 
2.032 d, for t/0987P=OGiP05209 2.094 d, and for 
F0231C=OGLP01092 1.802 d. 

We also analysed several merged data sets, containing 
about 90 - 100 Wh and F-corrected data points each. 
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but we were not able to improve the accuracy of the hght 
curve Fourier parameters of shorter P stars. In particular 
we cannot confirm the possible indications for the reso- 
nance effects in lOT pulsators reported in Paper II. The 
conclusion is that the accurate study of the light curves of 
Cepheids fainter than V ^ 22 requires a larger telescope 
than that adopted by us or by the OGLE project. 

Finally we mention the deep HST imaging of a field in 
the halo of IC 1613 performed by Dolphin et al. (|200lD; 



they found 11 short P Cepheids and 13 RR Lyrae stars. 



3. Period distribution 

The P distribution is shown in Figs. ^ ^ and ^ and is com- 
pared with those of other galaxies in the Local Group. In 
Fig. H the histogram of all the classical Cepheids is re- 
ported, while in Fig. || the stars are separated according 
to the pulsation mode; finally, Fig. || shows the homoge- 
neous distribution when the P of lOT pulsators is mul- 
tiplied by an average period ratio Pol Pi ^ 1/0.7. The 
adopted bin size is AlogP = 0.05; for IC 1613 in Figs. | 
andU it is AlogP = 0.1. Since we are interested in point- 
ing out the new results with respect to previous studies, 
the histograms are focused on intermediate and shorter 
P Cepheids (P ^ 50 d); a comment on longer P stars is 
however reported as a conclusive note in the Discussion. 

IC 1613. The histogram shown in Fig. || includes all 
the known 209 Cepheids, while the sample in Figs. ^ and 
^ is based on ou r surv ey and it includes the stars found by 
Dolphin et al. ( ^001 ). The samples should be considered 
almost complete for P longer than about 2 - 2.5 d (logP ^ 
0.4). Probably there are no lOT Cepheids with P 3 d. 

Milky Way. For the histogram in Fig. ^ we consid- 
ered 504 st ars in the Galactic Cepheid Database (Fernie 
et al. |l995| ). For Fig. | we have used the results on the 348 
stars with good light curves analyzed by the Brera-Merate 
group (e.g. Poretti [1994 , Antonello & Morelli 1996, and 
references therein), which include 36 identified lOT pul- 
sators. We recall that it is possible to discriminate between 
lOT mode and F mode stars using light curve Fourier 
parameters only when the P is shorter than about 5.5 d, 
since above this value the light curves tend to be similar. 
Therefore no identified \OT mode stars exist with longer 
P. On the other hand, the identified lOT Cepheids whose 
distribution is shown in Fig. || should constitute a fairly 
complete sample of these stars with P ^ 5.5 d in the 
vicinity of the Sun. 

LMC and SMC. The data for stars in LMC and SMC 



were taken from OGLE data archive (Udalski et al. 1999a, 
U dalski et al. |l999b| , Udalski et al. [19994 Soszynski et 
al 200C| ); they include also the double mode Cepheids 



(DMC). For SMC we have 1343 F stars (24 of them are 
DMCs pulsating in F and lOT mode), and 890 \OT stars 
(71 DMCs pulsating in lOT and 20T mode). For the 
LMC we have 791 P stars (20 DMCs pulsating in F and 
lOT mode), and 569 lOT stars (61 DMCs pulsating in 



lOT and 20T mode). Owing to OGLE observational con- 
straints, the sample does not contain stars with P longer 
than about 40 d; however this is not an issue for the prob- 
lem we will discuss here. On the other hand, the sample 
is fairly complete as regards the shorter P stars, and is 
well representative of the Cepheid population in the sur- 
veyed fields of the LMC and SMC. The distributions of 
LMC Cepheids show a long tail at shorter P, given by 
\OT pulsators. 

M31. After an analy sis of DIRECT M31 data archive 
(Mochejska et al. 1999[ and references therein), we pre- 
ferred to adopt the previous results on 392 stars reported 



in the GCVS (Samus [1995[ ). The reason is that DIRECT 
sample appears incomplete for P shorter than about 6 d, 
where the old data distribution shows a prominent max- 
imum; the difference is not due to possible different ob- 
served regions of the galaxy. As discussed by BIT, the 
sample produced by the photographic data should be con- 
sidered representative for logP ^ 0.7. We note in passing 
that there is not yet information available on the existence 
of lOT pulsators in this galaxy. 

M33. For M33 we adopted the results of the DIRECT 



s urvey published by Maori et al. (2001), Mochejska et al. 
( ^001a[) and Mochejska et al. (^001b[ ); the Cepheid lists 
were checked for the presence of the same star in differ- 
ent lists, on the basis of stellar coordinatae. The sample 
contains 654 stars. Given the adopted observational tech- 
niques, the Cepheids with P > 14 d are partially under- 
estimated. On the whole, keeping in mind this possible 
selection effect, we think the sample is representative of 
the Cepheid population for logP 0.5. Unfortunately the 
authors do not report the mode classification, even if un- 
certain; the PL relation for M33 shows the very probable 
presence of lOT mode pulsators with P up to 7 - 8 d. We 
have adopted the separation line reported in their Fig- 
ure 11 for discriminating between F and candidate lOT 
mode stars. Of course, the results of the analysis should 
be considered only as preliminary. 



4. Discussion and Conclusion 

We are trying to compare the frequency-period distribu- 
tion in different galaxies even if the samples have differ- 
ent degrees of completeness. The possible selection effects 
that hinder attempts to achieve complete samples were 
discussed by BIT, and we can revise them, taking into ac- 
count the progresses made in the last twenty years. Vari- 
ables with large amplitude variations arc more likely to 
be discovered than those with small amplitude variations. 
Since the amplitude is related to the P, the discovery of 
stars in the P ranges where the amplitude is larger is fa- 
vored. The photographic technique made it difficult to find 
Cepheids with amplitudes less than about 0.5 mag; for the 
CCD technique this limit is much lower (for example, this 
has allowed the detection of lOT pulsators in galaxies be- 
yond the Magellanic Clouds). The main selection effect 
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Fig. 2. Period distribution of Cepheids in galaxies of the 
Local Group 



is therefore the limiting magnitude of each survey, cou- 
pled with the detectable amplitude (which is related to 
the error of the measurements). This favors the selection 
of Cepheids with longer P since they are brighter, while 
the detectable Cepheids with shorter P and for a given 
pulsation mode tend to be those with larger amplitudes. 
In Sect. 3 we have reported for each galaxy the estimated 
shorter P limit, and a reliable comparison between galax- 
ies can be done above such a value. Finally, we mention 
the ever-present problem of discovering variables with P 
close to 1 d; this difficulty is related to the constraint of 
observing at one-day intervals. 




Fig. 3. Period distribution of Cepheids pulsating in F 
(dotted line) and lOT (continuous line) mode 



Many of the differences among the various galaxies in 
Figs. H, H and ^ can be interpreted in terms of differ- 
ent average metallicity. According to the data collected 



by van den Bergh (2000), the metallicity indicator values 



(for young population I stars) can be summarized as in Ta- 
ble 1. As is well known (see e.g. BIT), the location of the 
maximum of the frequency-period distribution increases 
with increasing metallicity. 

As regards IC 1613, taking into account the incom- 
pleteness of the sample, the distributions are in part sim- 
ilar to those of SMC. The maximum of the distribu- 
tion in Fig. ^ would be probably located at P ^ 1.6 d 
(logP ^ 0.2), and we should expect a larger number of 
stars with P ^ 2 d yet to be discovered. The peak at 
logP = 0.4 is intriguing, but we cannot give it much 
weight. It is interesting to compare the LMC and SMC 
m Figs. I and g after roughly correcting for the different 
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Fig. 4. Homogeneous period distribution for F and lOT 
mode Cepheids; the P of lOT pulsators was multiplied 
by a factor 1/0.7, which is an average value of the ratio 
Po/Pi (see text) 

Table 1. Values of metallicity indicators of galaxies of the 
Local Group (young stars; van den Bergh, 2000) 



galaxy 


[Fe/H] 


12 4- log (O/H) 


M31 




9.0 


Milky Way 


+0.06 


8.7 


M33 




8.4 


LMC 


-0.30 


8.37 


SMC 


-0.73 


8.02 


IC 1613 


-1.3 


7.86 



pulsation modes, the shape difference of the distributions 
(Fig. H) decreases significantly (Fig. U). On the other hand, 
the similarity of the shape of the distributions of lOT and 
F pulsators in LMC (Fig. ^) appears rather unique among 
galaxies. 

In spite of the different metallicities, the location of 
primary and secondary maxima of M31, Milky Way and 
M33 are practically the same. We can see in Fig. || that 
they are located at about 5 d (logP ~ 0.7) and 13 d 



(logP ^ 1.1; more exactly, between 12.6 and 14.1 d), re- 
spectively. The recent observations of M33 made by the 
DIRECT project give us the opportunity to update the 
study on the possible origin of the dip located at about 8 
- 10 d (logP ~ 0.9 — 1.0). According to the interpretation 



of Buchler et al. (1997), in the Milky Way a deficiency of 
F mode pulsators in this P range is expected owing to 
the instability of the pulsation cycle; since such stars will 
actually pulsate in the lOT mode, the deficiency must 
be compensated by an excess of lOT pulsators in the P 
range of about 5.6 - 7 d. As noted in the previous Section 
3 (Milky Way), it is not possible to discriminate the pul- 
sation mode for P > 5.5 d using light curve parameters 
alone, therefore for the present we cannot say how many 
lOT pulsators are actually located between logP ~ 0.75 
and ~ 0.84 in the Milky Way. In M33, according to the 
PL diagram reported by Mochejska et al. ( 2001b| ), there 
are some probable lOT stars located in such P range (Fig. 
^) , but it is difficult to affirm that there exists an excess of 
lOT stars; on the other hand also the dip, which is located 
near 10 d, is not very pronounced. The exercise shown in 
Fig. ^ indicates that indeed the dip tends to be filled by the 
lOT pulsators, if we change their P to that of the F mode 
one, but in any case there is still a bimodal distribution, 
i.e. there is an excess of stars with log P ~ 1.1 which must 
be explained. Finally, it is instructive to look at M31 even 
if there is no information on the lOT pulsators. In this 
case, even assuming a certain number of lOT pulsators 
with P within 5.6 - 7 d which more or less fill the dip, 
they will never cancel the evident bimodal distribution. It 
is important to remark th at we do not exclude the effect 
proposed by Buchler et al. (1997). We just say that it does 
not look sufficient to explain the bimodal distribution. 

It is difficult t o verif y which of the two interpretations, 
BIT or Chiosi's ( 1989 ) one, would be more acceptable. 
Roughly speaking, M33 data confirm the expectation of 
both that the secondary component depends on the metal- 
licity of the galaxy. If we make the reasonable assumption 
that the samples are fairly complete for a sufficiently long 
P (e.g. above 6 - 7 d), and we consider normalized dis- 
tributions of M33, Milky Way and M31, the data indi- 
cate that this component increases montonically with the 
metallicity. Then the question is which of the two hypothe- 
ses predicts this relationship. The answer requires a more 
detailed study. 

The clarification of some of these problems requires 
more accurate photometric data of Cepheids in the Local 
Group galaxies. Moreover, the longer P lOT Cepheids in 
the Milky Way have yet to be identified; a useful technique 
for such an identification could be the analysis of the ra- 
dial velocity curves. The preliminary results obtained by 
Moskalik & Ogloza ( 200C| ) suggest indeed that the Fourier 
parameters of radial velocity curves and also the phase lag 
between light and radial velocity curves could be mode 
discriminators for P ^ 5 d. This is a necessary study 
before beginning a reliable discussion of the Cepheids in 
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the Milky Way. Unfortunately, we think there is a further 
dilemma: assuming the excess of lOT pulsators does ex- 
ist, should it be the expected signature of the pulsational 
stability or would it be just the bimodal distribution of 
the lOT pulsators? 

As a final remark we note that, with respect to BIT, 
there are few new observational results on very long {P > 
100 d) Cepheids. The HST Key Project on Hq determina- 
tion was a good opportunity for a statistics on a large sam- 
ple of galaxies, but unfortunately the observations were 
limited to periods shorter than 100 d. According to BIT, 
Cepheids with logP > 2.0 are probably massive stars 
burning carbon in their cores, and those authors infer that 
a whole galaxy such as M31 would have to be searched to 
find just one Cepheid in the core-carbon-burning phase. 
However, these stars are present in dwarf irregular galax- 
ies such as LMC, SMC, IC 1613 an d, as r ecently observed, 
also in NGC 6822 (Antonello et al. ^00l[ ). BIT interpreted 
the discrepancy as further support for their suggestion of a 
two-component birthrate function, the second component 
giving 5-50 times as many massive stars as indicated by 
the one-component model. However, the nonlinear pul- 
sation characteristics could also play a role in this case. 
Aikawa & Antonello ( 2000| ) found that the nonlinear pul- 
sation cycle for P > 100 d is stable only for low metallic- 
ity values, therefore the very long P Cepheids should be 
preferably detected in metal-poor galaxies such as dwarf 
irregulars. This would be a simple explanation of what is 
actually observed. 
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